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Abstract 

We present a preliminary measurement of the composition of inclusive semileptonic B meson 
decays using 9.4 fb" 1 of e + e~ data taken with the CLEO detector at the T(4S I ) resonance. In 
addition to measuring the charged lepton kinematics, the neutrino 4-vector is inferred using the 
hermiticity of the detector. We perform a maximum likelihood fit over the full three-dimensional 
differential decay distribution for the fractional contributions from the B — > XJlv processes with 
X c = D, D* , D**, and nonresonant X c , and the process B — > X u lv. From the fit results we 
extract (M\ - M 2 D ) = (0.456 ± 0.014 ± 0.045 ± 0.109) GeV 2 /c 4 with minimum lepton energy of 
1.0 GeV and (M\-~M 2 D ) = (0.293±0.012±0.033±0.048) GeV 2 /c 4 with minimum lepton energy of 
1.5 GeV. The uncertainties are from statistics, detector systematic effects, and model dependence, 
respectively. 



* Submitted to XXI International Symposium on Lepton and Photon Interactions at High Energies, August 
11-16, 2003, Fermi National Accelerator Laboratory, Batavia, IL. 
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I. INTRODUCTION 



Nonperturbative QCD physics connects measurements of B meson decay properties to 
the properties of underlying weak flavor changing currents. Heavy-quark effective theory 
(HQET) combined with the operator product expansion (OPE) provides a framework in 
which many inclusive B decay properties can be calculated 

Using HQET, the inclusive B decay matrix elements can be expanded in powers of 
A-qcd/Mb- For each order in the expansion new nonperturbative parameters arise. At 
order Aqcd/Mb, there is Aand at order Aq CD /M b , there are Ai and A2. The nonpertur- 
bative parameter A relates the b-quark mass to the B meson mass in the limit of infinite 
b-quark mass. The parameter A2 is directly related to the mass splitting between the B* 
and B mesons. Moments of kinematic variables can be used to determine the A and Ai 
parameters. These parameters are properties of the B meson and are not specific to the 
decay mode being studied. This means that moments of the lepton energy and the X c mass 
in B — > XJu decays can be related to moments of the photon energy spectrum in B — > X s 7 . 
Furthermore, measurements of these parameters can be used in the calculation of absolute 
decay rates, providing an improved relationship between the B — > X c lu branching fraction 
and \V cb \. 

This paper presents a measurement of the first moment of the square of the hadronic mass 
in B — > XJu decays as a function of the lepton energy cut. Because the background from 
secondary leptons becomes large at low lepton energy, the lowest lepton energy considered 
is 1.0 GeV. Moments of the square of the hadronic mass have previously been measured 
by CLEO with a lepton energy cut of 1.5 GeV |2j and by BaBar with a lepton energy 
cut of 0.9 GeV [3]. As the lepton energy cut is made more restrictive, the terms that are 
higher order in Aqcd/Mb become more important, making the expansion less reliable 
There are significant complications in reducing the minimum lepton energy to 1.0 GeV. The 
backgrounds from continuum e + e~ — > qq events, fake leptons and secondaries all become 
more important. The structure of the nonresonant B — > X c lv decays also becomes more 
important. 

In order to reconstruct the recoiling hadronic invariant mass, it is necessary to reconstruct 
both the charged lepton and neutrino kinematics. The neutrino is reconstructed using the 
approximate hermiticity of the CLEO II and CLEO II. V detectors f6[ and the well known 
initial state of the e + e~ system produced by the Cornell Electron Storage Ring (CESR) 0. 
The experimental resolution on the neutrino energy has a narrow core with a full width at 
half maximum height of approximately 120 MeV and a broad tail of over estimation of the 
neutrino energy which extends up to 1.5 GeV. The square of the hadronic recoil mass, 
can be calculated if the B momentum, ps, is known, 

M\ = M 2 B + q 2 - 2E beam (E e + E v ) + 2\p B \\pi + Pv\ cos b4 u, 

where B .e v is the angle between the momentum of the leptonic system, q = pe + Pu, and 
Pb- Since the B mesons are the daughters of an T(4S) produced at rest, the magnitude 
of the B momentum is known and small, however its direction is unmeasured. The last 
term in the M\ equation depends on the unknown B momentum direction, but it is small, 
and therefore neglected in this analysis. Because of the neglected term, the M\ resolution 
depends on |g|. The kinematic variables q 2 = (pe + p») 2 and the pseudo-helicity angle of 
the virtual W, cos Owe, can also be calculated. The pseudo-helicity angle, cos Owe, is defined 
as the cosine of the angle between the lepton momentum in the virtual-W frame and the 
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virtual-W momentum in the lab frame. 

From the sample of events with a charged lepton and reconstructed neutrino, the dif- 
ferential decay rate is measured as a function of the reconstructed quantities q 2 , M|, and 
cos Owe. We fit the observed three-dimensional distribution to a sum of Monte Carlo models 
for the different hadronic final states: D, D*, D**, X c nonresonant, and X u . The models 
are constructed with Monte Carlo events generated with a full detector simulation [?J and 
have the same reconstruction cuts applied as the data. In addition to the B — > Xlv pro- 
cesses, there are backgrounds from fake leptons, continuum leptons and secondary leptons. 
The q 2 and cos Owi variables add information which contributes to the separation of the 
signal from the backgrounds, and the separation of the D and D* final states. The moment 
(M| — M D ) is calculated from the fit results for the final states containing charm using the 
same theoretical models as are used in the fit, but without any detector simulation. 

II. DATA SAMPLE AND EVENT SELECTION 

The data used in this analysis were taken with two configurations of the CLEO detector, 
CLEO II and CLEO II. V. An integrated luminosity of 9.4 fbT 1 was accumulated on the 
T(4S f ) resonance, E cm »s 10.58 GeV and an additional 4.5 fb _1 was accumulated at 60 MeV 
below the T(4S) resonance, where there is no BB production. Both detector configurations 
covered 95% of the 47r solid angle with drift chambers and a cesium iodide calorimeter. 
Particle identification was provided by muon chambers with measurements made at 3, 5, and 
7 hadronic interaction lengths, a time of flight system and specific ionization (dE/dx) from 
the drift chamber. In the CLEO II configuration, there were three concentric drift chambers 
filled with a mixture of argon and ethane. In the CLEO II. V detector, the innermost tracking 
chamber was replaced with a three-layer silicon detector and the main drift chamber gas 
was changed to a mixture of helium and propane. 

Events are selected to have an identified electron or muon with momentum greater than 
1 GeV/c and a well reconstructed neutrino. Additional criteria are used to suppress back- 
ground events from the e + e _ — > qq continuum under the T(45*) resonance. 

The identified leptons are required to fall within the barrel region of the detector (| cos#| < 
0.71, where is the angle between the lepton momentum and the beam axis). Electrons 
are identified with a likelihood-based discriminator which combines dE/dx, time of flight, 
and the ratio of the energy deposited in the calorimeter to the momentum of the associated 
charged track (E/p). Muons are identified by their penetration into the muon chambers. 
For momenta between 1.0 and 1.5 GeV/c, muon candidates are required to penetrate at 
least 3 interaction lengths and above 1.5 GeV/c, candidates are required to penetrate at 
least 5 interaction lengths. The absolute lepton identification efficiencies are calculated 
by embedding raw data from reconstructed leptons in radiative QED events into hadronic 
events. The rate at which pions and kaons fake leptons is measured by reconstructing 
Kg — ■> 7r + 7r~, D° — > K~ir + , and D° — *> K + tt~ using only kinematics and then checking the 
daughter particle lepton identification information. 

Neutrinos are reconstructed by subtracting the sum of the four-momenta of all observed 
tracks and showers not associated with tracks, V^bservedi fr° m the four-momentum of the 
e + e~ initial state, P^+ e ~, which is nearly at rest in the laboratory, 

Pu P e + e - V 'observed' 
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The errors made in this assumption are due to particles lost through inefficiency or limited 
acceptance, fake tracks and showers, and other undetected particles such as mesons, 
neutrons, or additional neutrinos. Several requirements are made to select events in which 
these effects are reduced and the neutrino four-momentum resolution is correspondingly 
enhanced. 

Because extra neutrinos are correlated with extra leptons, events with an identified lepton 
in addition to the signal lepton are rejected. The primary source of fake tracks is from 
charged particles which do not have sufficient transverse momentum to reach the calorimeter 
and therefore curl in the tracking chambers, returning to the beam axis. The portions of the 
track after the initial outbound portion may be accidentally included as additional tracks. 
Criteria have been developed to identify such errors and make a best estimate of the actual 
charged particles in the event. Events for which the total charge of the tracks is not zero are 
removed, reducing the effect of lost or fake tracks. Showers in the calorimeter associated with 
tracks in the drift chamber are not used, so as not to double count their energy. A neural 
network algorithm has been developed to reject secondary hadronic showers associated with 
showers that are already associated with tracks. 

A final neutrino reconstruction quality requirement is that the mass of the reconstructed 
neutrino to be small. The ratio of the reconstructed neutrino invariant mass squared, M^ iss , 
to twice the reconstructed neutrino energy is required to satisfy \M^ liss /2E u \ < 0.35 GeV/c 4 . 
This quantity is proportional to the energy of a lost or fake particle. After this cut, the 
reconstructed neutrino's energy is assigned to be the magnitude of the missing momentum, 
because the momentum is not dependent on the particle identification of the tracks and so 
has a better resolution than the direct energy measurement. 

Continuum events are suppressed by a combination of event shape and orientation criteria 
which exploit the fact that continuum events tend to be jet-like with a jet axis aligned with 
the beam axis, whereas BB events are more spherical and their orientation is uniformly 
distributed in the detector. The second Fox- Wolfram moment R2, of the energy flow in 
the event is required to be less than 0.4. In addition, a neural network is used to combine 
R2, the angle between the lepton and the event thrust axis, the angle between the lepton 
momentum and the beam axis, and the fraction of the total energy lying in nine separate 
cones around the lepton direction, which cover the full the Atc solid angle. The R2 cut is more 
than 99% and 95% efficient for B — > XJv and B — > XJv, respectively, while removing 60% 
of the continuum events. The neural net cut removes an additional 73% of the continuum 
background, while keeping 92% and 94% of the B — > XJv and B — > X u lv, respectively. 

After all cuts we observe 41411 events from CLEO II and 80440 events from CLEO II. V. 
The overall efficiency varies from 1.5% for B — > XJv nonresonant to 4.2% for B — > XJv. 

III. FIT TECHNIQUE 

The full three-dimensional differential decay rate distribution as a function of the re- 
constructed quantities q 2 , M\, and cos 9we is fit for the contributions from semileptonic B 
decay and backgrounds. The q 2 variable is replaced by q 2 j{E(_ + E v ) 2 for fitting purposes. 
The B decay modes are B — > Dlv, B — > D*lv, B — > D**lv, B — > XJv nonresonant, and 
B — > XJv. The backgrounds are classified as secondary leptons, continuum leptons, or fake 
leptons. A secondary lepton is a real lepton in a BB event whose parent is not a B meson. 
A continuum lepton is a real lepton in a continuum e + e~ — > qq event. A fake lepton is a 
non-leptonic track from either a BB or a continuum event which is identified as a lepton. 
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The normalization of the continuum lepton component is determined from the data taken 
below BB threshold. The normalization of the fake leptons is determined from the measured 
fake rates and the measured track spectra. The contributions of these two backgrounds are 
therefore not allowed to vary in the fit, while those of the secondary leptons and all of the 
B — > Xlv modes are. 

We perform a binned maximum-likelihood fit where component histograms are con- 
structed from weighted Monte Carlo or data events. The fit uses electrons and muons 
simultaneously, with a separate set of histograms for each. The likelihood is implemented 
to take into account the histogram statistics using the method described in reference (§[. 
Projections of the Monte Carlo simulations of reconstructed quantities q 2 , M\, and cos 6^ 
for the various B — > Xlv modes are shown in Figure H Projections of the data and fit result 
are shown in Figure 121 




FIG. 1: Monte Carlo simulation of the distributions of the reconstructed quantities q 2 , M^, and 
cos9\yi for the various B — > Xlv modes. The modes are B — > X u lv (solid), B — > Dlv (short dash), 
B — > D*lv (dots), B — > D**lv (dot-dash), and B — > X c lv nonresonant (long dash). The curves 
are normalized to have unit area to facilitate comparison of the shapes. Note that due to finite 
resolution My can be less than zero. 



The B — > Xlv modes, secondary leptons and real leptons from the continuum are modeled 
with events from a GEANT l7f simulation of the CLEO detector that are reconstructed in 



6 




^ o 12000 

> 

(51 0000 



H 8000 

c 

<u 
> 

W 6000 



4000 



2000 



'i 1 1 1 i 1 1 1 1 



• 




10 15 20 25 30 

q 2 GeV 2 



i ... i ... i 



-4 -2 2 4 6 8 10 12 14 



GeVV 




-0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 

Cos0 w , 



1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 

Lepton Momentum GeV/c 




0.5 



1 1.5 2 2.5 3 

Neutrino Energy 



3.5 4 

GeV 



□ b^u 

□ Dlv 
^ D*lv 

** 

Dlv 

□ Nonresonant-X c lv 

□ Backgrounds 



FIG. 2: 



and M|, 



Projections of the fit results in the reconstructed quantities, on the first row are q 2 
on the second are cos6y/n and E#, and on the third row is E v . From the top in each 
figure the histograms are the contributions from B — > X u lu (barely visible shaded), B — > Dlv 
(open), B — > D*lu (sparsely hatched), B — > D**lv (densely hatched), B — > X c lu nonresonant 
(light shaded), and the sum of the backgrounds: continuum, secondary leptons, fake leptons (dark 
shaded). The step in the lepton momentum distribution is due to the looser muon identification 
used below 1.5 GeV/c for which there is a higher fake rate and higher efficiency. 
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the same manner as data events. The B — > Dlv and B — > D*lv modes are simulated with 
HQET using the PDG 0] averages of measurements of the form factors. The B — > D**lv 



and 5 — * X^iz/ modes are simulated using form factors from the ISGW2 model The 



X c nonresonant modes are simulated with the Goity and Roberts model l3| in which the 
D and D* contributions are excluded. 

The fake leptons are modeled with data events where a track is selected to be treated as a 
lepton. The events are then unfolded bin by bin to extract the 7r and K contributions, which 
are then multiplied by the measured fake rates. This models fake leptons from both BB 
and continuum. This method also provides an absolute normalization for the fake lepton 
contribution to the data sample. The real leptons from continuum are modeled with a Monte 
Carlo simulation which has been tuned to replicate the appropriate charm spectra; charm is 
the source of most leptons from continuum. The models of both continuum and fakes have 
been validated and constrained by comparisons with the 4.5 fb _1 of off-resonance data. The 
secondary leptons are modeled with CLEO's generic BB Monte Carlo which has also been 
tuned to replicate measured charm spectra and semileptonic charm decay measurements. To 
address the issue of the poorly understood charm counting in B decays, the charm counting 
in the BB Monte Carlo simulation is tuned to saturate the theoretically predicted level of 
charm production 

The B — > Xlv models were varied to assess the dependence of the results on the model. 
The B — > Dlv and B — > D*lv modes were varied within the range of the errors on the 
form factor measurements. The curvature of the form factors has not been measured and 



is usually constrained by theoretical predictions [14j. Because the data have an excess 
above the model in the q 2 region between 5.0 and 8.0 GeV 2 /c 4 , the curvature was set to 
50% of its predicted value and varied by ±50% of the prediction. The B — > D**lv form 
factor was replaced by a model inspired by HQET calculations (l5j |. The slope of the 
B — > D**lv and B —>■ XJv nonresonant form factors in the q 2 dimension was also varied. 
The dominant model dependence is from the hadronic mass distribution of the B — > XJv 
nonresonant mode. This is conservatively reweighted with a series of Gaussians restricted 
to the kinematically allowed region. The means of the Gaussians are allowed to range from 
M D + M w to 3.5 GeV/c 2 with variances ranging from 0.25 GeV 2 /c 4 to 1.25 GeV 2 /c 4 . The 
B — > X u lv simulation is varied from an all-nonresonant model to the nominal ISGW2 model 
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with a hybrid of the two in between. The all-nonresonant model differential decay 
rate corresponds to the prediction of HQET combined with CLEO's measurement of the 
B — > X s 7 spectral function r| 17 1. The maximum deviation of the B — > XJv nonresonant 
mass Gaussians is added in quadrature with the deviation of the other model variations to 
get the total model dependence. 

Radiative corrections play an important role in the measurement of the M\ distribution. 
The reconstructed M\ is defined to be the mass squared of the system recoiling against the 
charged lepton and the neutrino. If a photon is radiated by the lepton in the event, it will 
be included in this definition of the recoil system, 

Mi = (^-^-^) 2 = K+^) 2 - (1) 
The goal of this analysis is to measure the mass squared moment of the recoiling hadronic 
system p 2 , not its combination with the radiated photon (ph+P-y) 2 '• To correct for this effect, 
the data are fit using fully simulated GEANT Monte Carlo events in which the PHOTOS 
package 0] has been used to generate radiated photons. The moments are calculated from 
the fit results and the theoretical models of the hadronic mass distributions for each mode 
and thus do not have a shift due to the radiative corrections. 
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However, knowledge of the radiative corrections is not complete. PHOTOS implements 
an algorithm based on a splitting function which applies the same physics at 0{a) as the 
prescription of Atwood and Marciano 19] . PHOTOS also includes a prescription for forcing 
the kinematics of the decay to conserve momentum in addition to energy. The implementa- 
tion of PHOTOS has been checked in detail against a private implementation of the same 
algorithm. The applicability of the PHOTOS algorithm is however not exact. The calcu- 
lation ignores the internal structure of the hadronic system. Richter-Was |2(| has made 
a comparison of the PHOTOS and Atwood-Marciano prescriptions with an exact order-a 
calculation of the radiative corrections to the B~ — > D°£~V differential decay rate and has 
found agreement at the 20% and 30% level, respectively. Because we must extrapolate to 
the other B — > XJu modes, we make a conservative estimate that the PHOTOS calculation 
can be trusted to ±50%. 

The application of the radiative corrections is further complicated by the fact that the 
generated photons are low energy and often lost. When the photon is lost it can cause 
the event to fail the missing mass cut. If the event does pass the missing mass cut, the 
reconstructed neutrino will be biased toward high energy, pushing the reconstructed M\ 
toward the true hadronic mass squared without the photon. If neglected, this would increase 
the measured (M% — M D ) moment with a 1.0 GeV lepton energy cut by 0.082 GeV 2 /c 4 , 
before detector effects are included. This would be reduced to 0.037 GeV 2 /c 4 after detector 
effects. 

The method of neutrino reconstruction adds a large amount of kinematic information to 
each event. However, it also adds significant potential for systematic errors. The resolution 
on the neutrino kinematics is affected by the models of the signal, the other B in the 
event, and the detector response. The GEANT Monte Carlo simulation does not perfectly 
reproduce the track and shower efficiencies and fake rates, nor are B decays well enough 
understood that the inclusive particle distributions are well known. For this analysis we 
employ a reweighting method in order to quantify the effects of these uncertainties on our 
results. For example, to study the effect of the tracking efficiency uncertainty, the Monte 
Carlo events in which tracks are lost are given a higher or lower weight in constructing the 
component histograms for the fit. The scale of the variation can in general be constrained 
by direct measurements of the quantity being varied. One important example is that having 
a K\ in the event adds a tail to the neutrino resolution. The inclusive Kg spectrum in BB 
events has been measured and can be used to constrain the K\ spectrum. The dominant 
detector systematics are fake showers (generally splitoffs from hadronic showers), final state 
radiation, and fake leptons. The detector systematics are summed in quadrature to arrive 
at the overall detector systematic error. 



IV. RESULTS 

The fit results can be used to calculate a branching fraction for each of the hadronic final 
states B — > Dlu, B — > D*lu, B — > D**lu, and B — > XJu nonresonant. From these branching 
fractions the moment (M^ — M D ) with a certain lepton energy cut can be calculated as 
follows: 

(M* x -MlU= ^™ mmC f m , (2) 
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where m m is the moment of mode m with the cut applied, c m is the fraction of mode m 
passing the cut, and B m is the branching fraction of mode m. The quantities m m and c m 
depend only on the model. The measured branching fractions, B m , depend on the model, 
the detector simulation, and the data. 

For each lepton energy cut, the fit is repeated with the cut applied using the reconstructed 
(lab frame) quantity. Theoretical calculations of (M\ — M D ), however, are made with the 
cut in the B rest frame 0, 0] . The moments reported are therefore calculated in the B rest 
frame. This involves a small extrapolation. 

The resulting moments are shown in Table |U These moments are highly correlated 
because they share some fraction of the data and also use the same models and detector 
simulation. The statistical correlations are assessed with a simple Monte Carlo simulation. 
The model dependence and detector systematic correlations are assessed by summing the 
covariance matrices from each systematic error. The resulting covariance matrix is shown 
in Table |HI An alternative representation of the data is to use the {M% — M D ) moment 

with a 1.5 GeV lepton energy cut and the difference between the {M\ — M D ) moments at 
1.0 and 1.5 GeV. These two measurements are significantly less correlated. We find 

(Mj^^oGeV - (M|) B<>1 .5GeV = (0.163 ± 0.014 ± 0.036 ± 0.064) GeV 2 /c 4 

and a covariance of this value with the (M|— M D ) E t >i.5GeV 

moment of 2.242 x 10~ 3 GeV 4 /c 8 . 

TABLE I: {M^ — M D ) versus the lepton energy cut. The errors on the entries in the table are the 
statistical, detector systematics, and model dependence, respectively. 



Cut (GeV) 


(Ml-M 2 D ) (GeV 2 /c 4 ) 


Ee 


> 1.0 


0.456 ± 0.014 ± 0.045 ± 0.109 


Ei 


> 1.1 


0.422 ± 0.014 ± 0.031 ± 0.084 


Et 


> 1.2 


0.393 ± 0.013 ± 0.027 ± 0.069 


Ei 


> 1.3 


0.364 ± 0.013 ± 0.030 ± 0.054 


Ei 


> 1.4 


0.332 ± 0.012 ± 0.027 ± 0.055 


Ei 


> 1.5 


0.293 ± 0.012 ± 0.033 ± 0.048 



TABLE II: Covariance of (M^—M D ) moments with different lepton energy cut. The rows/columns 
correspond to progressively more restrictive lepton energy cut starting with 1.0 GeV and ending 
with 1.5 GeV. 

/1.421 1.096 0.913 0.724 

1.096 0.864 0.719 0.577 

0.913 0.719 0.612 0.495 

0.724 0.577 0.495 0.413 

0.735 0.574 0.499 0.418 
\0.618 0.486 0.422 0.352 

Figure H3 shows a comparison of the results with previous measurements [2j, |3( and a 

2 

prediction based on HQET constrained by this measurement of (M\ — M D ) E t >i.5 and the 



0.735 0.618\ 

0.574 0.486 

0.499 0.422 

0.418 0.352 

0.447 0.373 

0.373 0.324/ 



x 10^ 2 GeV 4 /c 8 
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first B — > X s 7 photon energy moment The results agree well with both. The theory 
bands shown in the figure reflect the experimental errors on the two constraints, the variation 
of the third-order HQET parameters by the scale (0.5 GeV) 3 , and variation of the size of 



the higher order QCD radiative corrections [21 
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FIG. 3: The results of this analysis compared to previous measurements 0, EJ and the HQET 
prediction. The theory bands shown in the figure reflect the variation of the experimental errors 
on the two constraints, the variation of the third-order HQET parameters by the scale (0.5 GeV) 3 , 
and variation of the size of the higher order QCD radiative corrections 21 ] . 



V. CONCLUSION 



We have measured the moment {M\ — M D ) as a function of the lepton energy cut with 
a minimum cut of 1.0 GeV. We obtain the preliminary result, 

(M 2 X - M 2 D )E e >i.oGeV = (0.456 ± 0.014 ± 0.045 ± 0.109) GeV 2 /c 4 . 

This result provides a test of the HQET predictions using a larger portion of the available 
decay rate. We have also improved the precision of the measurement with a 1.5 GeV lepton 
energy cut, obtaining the preliminary result 

(M| - M 2 D ) Ee>1 . 5GeV = (0.293 ± 0.012 ± 0.033 ± 0.048) GeV 2 /c 4 . 

Measurements such as these can be used to improve the determination of the Standard 
Model parameter \V c b\- They also provide a valuable test of the HQET inclusive predictions 
which are used in the extraction of \V u b\ and \V c b\ and in the interpretation B — > Xgj . 
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